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Abstract 

We argue that the minimal chiral background for the two-pion exchange 
nucleon-nucleon interaction has nowadays a rather firm conceptual basis, 
which entitles it to become a standard ingredient of any modern potential. In 
order to facilitate applications, we present a parametrized version of a con- 
figuration space potential derived previously. We then use it to assess the 
phenomenological contents of some existing NN potentials. 
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I. INTRODUCTION 



Nowadays there are several nucleon-nucleon potentials which reproduce experimental 
information with high precision (l|-|6| . In all of them, the long and medium range interactions 
are associated with the exchanges of one and two pions, whereas there is a wide variation in 
the way short distance effects are treated. This region is theoretically uncertain and good 
fits to the data usually require many free parameters. The availability of alternative models, 
although welcome, points out to the need of criteria for discriminating among the various 
potentials. 

All potentials become progressively less reliable as one goes from large to short dis- 
tances. The pionic tail of the interaction is quite well established and its only source of 
uncertainty is the value of the pion-nucleon coupling constant. It determines many observ- 
ables of the deuteron and large angular momentum scattering |7|-[T0|]. As far as the next 
layer is concerned, one finds some consensus in the literature regarding the importance of 
delta intermediate states and pion-pion correlations. Nevertheless, the treatment of details 
of these interactions is not uniform, a fact reflected in the existence of different profile func- 
tions for a given component of the potential. It is in this framework that chiral symmetry 
may prove to be useful and provide guide-lines for a proper assessment of the various existing 
potentials. 

Chiral symmetry corresponds to the idea that strong interactions are approximately 
invariant under transformations of the group £77(2) x SU(2). This symmetry was developed 
in the sixties, in the framework of hadron physics, but the reasons behind its successes 
became clear only after the formulation of QCD, the theory decribing the basic interactions 
of quarks and gluons. In QCD the quark masses of the low-lying £77(2) multiplet are 
very small and the Lagrangian becomes chiral invariant when they are neglected. QCD 
predictions can be directly tested only at high energies, where perturbation makes sense. At 
low energies, on the other hand, its non-Abelian structure makes calculations very difficult, 
and one is forced to resort to effective hadronic theories. In order to remain as close as 
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possible to the fundamental level, the hadronic Lagrangian must have the same symmetries 
as QCD, even the approximate ones, such as chiral symmetry. 

One of the major successes of chiral symmetry at hadron level were the predictions 
made for low-energy 7rN scattering. Early chiral calculations of this process considered just 
minimal systems, containing only pions and nucleons [|TT],[nj. Motivated by phenomenology, 
this model was extended and one has learned that chiral amplitudes at tree level, including 
nucleon and delta poles, rho exchanges and corrections associated with 7rN sigma term do 
reproduce satisfactorily subthreshold and scattering data up to 300 MeV. 

Chiral symmetry at hadron level may be implemented by means of either linear or non- 
linear Lagrangians. The fact that no serious candidate has been found for the a meson 
favours the latter type of approach. There are two forms of non-linear Lagrangians which 
are especially suited for the 7rN system. One of them is based on a pseudoscalar (PS) 
7rN coupling supplemented by a scalar (S) interaction, equivalent to the exchange of an 
infinitely massive a meson, and denoted as PS + S scheme (Eq. [I]). The other one employs a 
pseudovector (PV) 7rN interaction and a vector (V) term, which could represent the exchange 
of an infinitely massive p meson, constituting the PV + V scheme (Eq. |2|) . 

£ps+s = # (J II ~ <P 2 + ir ■ 07 5 ) ip , (1) 

1 — q — 

4/; 2m 
Both approaches yield the very same amplitude for the 7rN scattering when the axial coupling 

constant Qa is equal to 1. The extention to the case ^ 1 can be done quite naturally 

in the PV + V approach, and a little less so in the PS + S case. The fact that physical 

results should be independent of the representation used to implement chiral symmetry was 

discussed in very general terms by Coleman, Wess, and Zumino [L3||. In the case of 7rN 



interactions, this point was emphasized by Coon and Friar [II . 

Chiral symmetry has very little to say about the one pion exchange NN potential since, 
in the PS + S approach, predictions for this component do not depend on the scalar inter- 
action. Moreover, PS and PV couplings are equivalent. The intermediate range part of the 



interaction, on the other hand, is due to the exchange of two pions and closely related to 
an off-shell 7rN amplitude. Thus, as pointed out by Brown and Durso in 1971 ||15|| , chiral 
symmetry is quite relevant in this region. In the first calculations of this component of the 
potential, chiral symmetry was included disguised in the form of soft-pion theorems fT5| , |l6|| . 
Nowadays one has realized that the use of effective field theories makes the implementation 
of chiral symmetry into hadronic systems much easier. 

In the early seventies, approaches to the intermediate range NN interactions based on field 
theory disconsidered chiral symmetry. Nevertheless, they were very important in setting the 
method to tackle the problem. In 1970, Partovi and Lomon [lTj performed a paradigmatic 



calculation based on Feynman diagrams, using just PS couplings without chiral countert- 
erms. A little later, Partovi and Lomon [|18j included the a meson in their calculation, used 



chiral symmetry to determine the crN coupling constant, but did not group diagrams into 
chiral families, producing a potential with an uneven number of loops. About ten years 
later a rather comprehensive study was performed by Zuilhof and Tjon |I9fl , who compared 
direct and crossed box diagrams obtained using both PS and PV couplings, without paying 
attention to chiral symmetry. 

The first modern field theoretical calculation of the two-pion exchange nucleon-nucleon 
potential (7T7tE-NNP) with strong emphasis on chiral symmetry was produced by Ordonez 
and van Kolck p0|. They applied the conceptual pertubative framework developed by 



Weinberg [^TJ in a systematic discussion of corrections to pure OPE interactions, associated 
with form factor effects, exchange of two pions, and three body forces. The problem was 
formulated using PV coupling for the pion, and expressions were given in momentum space. 



In 1992, Celenza, Pantziris and Shakin [22] studied the influence of chiral symmetry 



on the NN interaction at one-loop order. These authors considered both PS and PV pion- 
nucleon couplings in order to derive the qualitative features of the scalar-isoscalar component 
of the interaction. They employed tensor decompositions of diagrams in momentum space, 
and one of their main conclusions was that the irreducible isoscalar interaction is attractive. 
In a series of independent works, other authors also addressed to the same problem, 



namely the determination of the role of chiral symmetry in the two-pion exchange NN 
interaction. A momentum space evaluation of the scalar-isoscalar component was produced 
by Birse |23[], in the framework of the linear a-model. Friar and Coon |[24fl , used a more 



general Lagrangian in order to obtain a momentum space potential for infinitely massive 
nucleons. Finally, in another calculation, we used non-linear Lagrangians based on 
covariant derivatives and obtained a non-relativistic potential in coordinate space. 

All these calculations fUJH^II are based on chiral symmetry and should, in spite of 



technical differences, yield the same predictions for the two-pion exchange potential. How- 
ever, inspection of the conclusions presented by the various authors shows that this is not 
the case. Discrepancies may be traced back to two different sources, namely the methods 
employed for regularizing divergent integrals and the procedures used to subtract the iter- 
ated OPEP. We here consider only the latter problem. It is well known that the meaning 
of an iteration is defined by the dynamical equation adopted in a given problem. Thus, the 
use of the relativistic Bethe-Salpeter equation [^] or the non-relativistic reductions due to 
Brueckner and Watson |27| or Blankenbecler-Sugar J2S| correspond to different forms for the 
potential. An explicit discussion of the non-relativistic case can be found in Ref . , where 
it is shown that physics is independent of the particular approach employed. 

Another important feature of the problem is that physical results in a consistent calcu- 
lation should also not depend the particular representation used in the implementation of 
chiral symmetry [0,0 • In our calculation [^] , we have shown explicitly that both PS + S 



and PV + V Lagrangians correspond to the same non-relativistic potential. Of course, there 
is no reason to expect that things would be different in the case of the Bethe-Salpeter equa- 
tion. However, we believe that the explicit demonstration of this property in the relativistic 
framework may prove to be non-trivial. The basic idea of the Bethe-Salpeter equation is to 
decompose the full amplitude into two-particle propagators and irreducible kernels. In the 
PS + 5* case, this kind of decomposition poses no special problem. The PV + V case, on 
the other hand, is tricky, since its box diagram contains contributions in which the pole of 
the nucleon propagator is canceled and hence must be interpreted with great care. Cancel- 



lations of nucleoli poles are, in fact, already present in the intermediate 7rN amplitude that 
determines the 7T7tE-NNP p5||29[j . Therefore, in the present state of the art, we believe that 
non-relativistic potentials provide a much safer ground for the study of chiral symmetry. 

The minimal chiral model for the 7T7tE-NNP, involving only pions and nucleons, has 
nowadays very solid conceptual foundations. It has attained the same status as the OPEP 
in the late sixties and one may expect it to become a necessary ingredient of any modern 
NN potential. On the other hand, this minimal chiral model fails to reproduce experimental 
information in the case of the intermediate 7rN amplitude incorporated into the NN inter- 
action. Thus, in the 7T7tE-NNP, the minimal model corresponds to a contribution which is 
necessary, but not sufficient to describe reality. This requirement demands the inclusion of 
other fields and processes in the model. An extension, in the framework of chiral symmetry, 
has been recently considered by Ordonez, Ray and van Kolck 3(J, with promissing results. 

In all calculations of the chiral background to the 7T7tE-NNP, results are given in form of 
cumbersome integrals, which must be evaluated numerically. In order to facilitate applica- 
tions, in Sec. II of this work we present a parametrized version of our configuration space 
potential, valid up to 0.5 fm and completely free of cut off parameters. The knowledge of 
the minimal chiral 7T7tE-NNP may be used to assess the phenomenological content of a given 
potential. This information can be obtained by subtracting the OPEP from it and then 
comparing the remainder with the chiral two pion exchange background. In Sec. Ill, we 
consider the dTRS 0, Paris ||, and Argonne || potentials in order to produce an instance 
of this kind of study. 



II. PARAMETRIZATION 

Our calculation of the 7T7tE-NNP is based on the Blankenbecler-Sugar reduction of the 
Bethe-Salpeter equation. In the PS + S scheme, its dynamical content is associated with 
the five diagrams displayed in Fig. [I]. Therefore we label the corresponding individual 
contributions by Q , A , N and □, where the last one also includes the subtraction of the 
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iterated OPEP. It has the general formQ 



V(r) = [(Vf + Vg) + LS (V LS + V™)} + 

(3 + 2r« ■ r^) [v£ + O ss V^ + 6 LS V£ S + 6 T V^} + 
(3 - 2t« ■ t®) [VS + OssV* 3 + LS V n LS + 6 T v£ 



(3) 



where the spin operators are given by Oss = cr ^ 1 ' ) • cr ^ 2 \ Ols = L ■ \{cr^> + cr^ 2 )}, and 
Ot = 3cr( 1 - ) • f cr^ 2 ) • f — <j^> • cr^ 2 \ whereas (r™ and represent spin and isospin matrices 



for nucleon (z). 

In the case of the bubble diagram, the leading contribution to the asymptotic potential 
can be calculated analytically, as shown in appendix. This results sets the pattern for the 
parametrization of the other components of the force. 

Our numerical expressions represent the various components of the potential in MeV, 
and are given in terms of the adimensional variable x = fir, where fi is the pion mass. We 
keep the 7rN coupling constant g as a free parameter and adopt the values \x = 137.29 MeV 
and m = 938.92 MeV for the pion and nucleon masses respectively. 

In general, the parametrized expressions reproduce quite well the numerical results of 
Ref. [KJ, except for a few cases and regions where the discrepancies become of the order of 
0.25%. Our results are listed below. 



The profile function for the central potential has the following common multiplicative 
expression 



A. Central Potential 




(4) 



The parametrization of each diagram gives: 



■ee Eq.(59) of Ref. @. 
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(~* , v n/ J 51.0923 6.54068 1.26190 0.1307061 
V C x = F c (x) -275.364 + — + — 5 

u I X X 2 X 6 X 4 J 

Vg(x) I'A.r) (3 13.558 ll ° J 



+ (135.249 + 14.6514 • x + 6.43825 • x 2 ) ■ e" - 397835 ^} (6) 

YL \ \/ X X 



577.210 378.004 133.374 19.5061] 

+ — 3 Jc- + Is (7) 



Xv^X X 2 X 2 y / X X 3 , 

V?(*) = F C (X) (-25.9987 + - 2^Z^1\ . e -[0.101214., + 0.00123687^] (g) 

I X X 2 J 



B. Spin-spin potential 

The multiplicative factor for the spin-spin potential is the same as that of the central 
potential, and receives contributions from the box and crossed diagrams only, which are 
given by 

1.05042 0.421043 0.0284309 



r^/xf 1.05042 0.421043 0.0284309 nMrnnn ,™ 4t 1 
x = FJx) \ 0.408084 + + 0.215829 • e" 1 - 2344 -* L 9 

I X X 2 X 6 J 

V*>(x) = F„(«) {0.399845 + i°™ + _ °- 03 °f 71 + 0.0371333* . e" } (10) 

I XX 2 X 3 J 



C. Spin-orbit Potential 

The spin-orbit multiplicative function is 

'-M-(£)Vs)£^ «") 

and individual contributions are: 

5.51078 0.994157 0.217562 0.0336067 0.002446201 



TO/ >. . . r 5.51078 0.994157 0.217562 0.0336067 0.002446201 , . 
V LS (x) = F LS (x) {-5.88744 - — — + — — + } (12) 

V^ 5 (x) = F L5 (x) {7.34548 + 



2.15233 0.381025 



x x 2 
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-0.419984-rr 



+ (7.48798 - 0.448484 • x + 0.391431 • x 2 ) ■ e 
14.7389 12.6506 6.27374 1.66637 + 0.1842641 

rr* / ry rv>2i rf 2 / ,-y» rv> 3 i-y»3 / <-y> I 

XyX X X V/ X X X V ^ I 

„ r ,. N n ,,f „ 1.58089 0.319790 0.01964611 . 

V^ s (x) = F LS (x) -1.19527 + • 1- 0.0160259 • x] 

I X X 1 X 6 J 



-1 



D. Tensor Potential 



The commom factor for the tensor potential is 

F ^ x) = (£) 4 i 1 + 1 + h) U ■ (16) 

It receives contributions from the box and crossed diagrams only, which have the form 

rp f . n / J 0.510720 0.0597556 
V M T x = F T x -0.204041 + 

I x x z 



+ (1.32932 - 0.939553 ■ x + 0.706050 ■ x 2 ) ■ e" 2 - 29686 '*} 



0.521123 0.352463 0.135028 0.0303012 0.00297840' 



Va (x) = F T (x) -0.246349 + + 



nr* rf 2 ,y»3 i -y»4 /y»5 

.X' lX-' J, X 



III. PHENOMENOLOGY 

As discussed in the introduction, one may expect a realistic two-pion exchange nucleon- 
nucleon potential to be composed of a minimal chiral background supplemented by other 
terms determined by phenomenology. In this section we estimate crudely the relative weights 
of these two types of contributions, by comparing the parametrized expressions of the pre- 
ceding section with tree realistic potentials available in the litterature, known as dTRS 0, 
Paris @, and Argonne vl4 ||. 

In the minimal chiral potential the 7rN coupling constant may be kept as a free parameter, 
whereas in the realistic versions considered it assumes slightly different values. Therefore, in 
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order to compare properly the various results, we subtract the OPEP contribution from each 
potential and then divide the remainder by the fourth power of the 7rN coupling constant 
used in it. We then consider, for various channels, the ratio of this quantity by that obtained 
from the parametrization of Sec. II with g — 1. 

The ratios for the (T = 1, S = 0) central potential are displayed in Fig. |2|. For distances 
greater than 7.5 fm, the ratio of the Argonne potential becomes a constant since, as the 
minimal chiral background, it behaves asymptotically as Yukawa functions with double pion 
mass. In the case of the dTRS potential, the ratio vanishes because it is relatively short 
ranged. The tail of the ratio of the Paris potential, on the other hand, is rather large and 
not stable up to 10 fm, probably due to the way it is parametrized. 

The realistic potentials considered here reproduce well NN observables due to their be- 
haviour in the region below 2.5 fm. Inspecting Fig. Q one note that, in spite of important 
differences in this region, the three potentials display a sort of rough coherence. In all cases 
the fit of observables seem to require contributions with the same sign and about four times 
bigger than the chiral background. The case of (T = 1, S = 1) central component is similar. 
The rough coherence among the realistic potentials is also present in the (T = 0, S = 0) 
and (T = 0, S = 1) central channels, but the ratios in the region between 1 and 2.5 fm are 
negative. 

In Fig. |H we show the profile function for V c , the spin-isospin symmetric component of 
the central potential, which corresponds to the combination V c = + V£ + 3V£ + 3V^% 
and is relevant for symmetric nuclear matter. The chiral background, in the physically 
relevant region, is much less attractive than the realistic forces. 

In Fig. |] we display the ratio of the (T = 0, S — 1) spin-orbit potentials. In this 
case there is no coherence among the realistic potentials, even regarding their signs. Both 
the tails and medium distance behaviours of the various ratios differ significantly from the 
background. The (T — 1, S — 1) component behaves similarly. 

Finally, the ratios for the (T = 1, S = 1) tensor potential are given in Fig. ^ where is 
possible to note that the realistic potentials are roughly coherent below 2.5 fm and about 3 
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times greater than the background. In the case of the (T = 0, S = 1) component, on the 
other hand, this ratio is about 5. 

In summary, the crude comparisions made in this section allow one to conclude that 
the minimal chiral background for the two-pion exchange nuclear potential accounts for less 
than 25% of the values need in order to reproduce observations. 
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APPENDIX: ASYMPTOTIC BEHAVIOR OF THE BUBBLE DIAGRAM 

The asymptotic behaviour of the bubble diagram [Fig. [TJ(a)] can be calculated analyti- 
cally. Its contribution to the potential is given by |25| 

t 
4m 2 

where 



Vo=<6^J«./(»).J 1 , (Al) 



r d*k 1 
1 ~ J (2^) 4 (P- A i 2 )(A;' 2 - A i 2 ) ' ( ' 

and the operator ■ has the following non-relativistic spin structure 

JW ■ 1^ = 1 _ °™ , (A3) 
2m z 

Ols being the spin-orbit operator. 

After the usual integration, the function J\ may be written as 

Ji(A) = (A 2 - /x 2 ) f da dP-—- 1 , (A4) 

(4:7v) z Jo Jo A 2 a(l — a) + — p) + pA z 

where A is the momentum exchanged in the interaction and A is a monopole cut off param- 
eter. The corresponding expression in configuration space is given by 

A 2 — n 2 r 1 r 1 ^ a 1 e -Mir 

= 1 ia \s 4 / da / d Pr— m — ? — ™ ' ( A5 ) 

(4^)-^ Jo Jo (a + p) — (a — p) z r 
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where 



Mi = 2 



1 -a-/3)A 2 + (a + (3)fi 2 



1/2 



(A6) 



(a + (3) -(a- f3) 2 

In the limit A — > oo, the integrand is non- vanishing only in the neighbourhood of the 
point (1 — a — f3) — 0. That means that, except in the term proportional to A 2 , we may set 
/3 — 1 — a. With this approximation, one has 



A 2 — u 2 r l ri—a 



J_ 2[(l-a-^-^] 1/2 .i 



D 2 



where 



(A7) 



D = 1 - (2a - l) 2 . 



(A8) 



Defining a new variable s by 



s 2 = (1 - a - (3)A 2 + i? , 



(A9) 



one has 



g rl I ro 

JAr) = i- — — / da —: — / 
U ; 4tt 3 Jo D 2 rJ„ 



ds se o . 



(A10) 



(4tt) 3 Jo D 2 

In deriving this result we used the fact that A is very large. Performing the integration and 
using the variable x = fj,r, one obtains 

1 1 



4u 3 f 1 
Ji(r) = i — — I da 



(47r) 3 Jo 



Dx 2 2 x 3 



- 2£ 

e d . 



(AH) 



For very large values of x, the main contribution to the integral comes from the region 
around a — |, where D is maximum. In order to explore this property, we use the variable 
9 = 2^fx[a — |), expand D around 9 = and write 



JAx) 



2/x 3 



-2x 



d9 



1 + J_ f 1 + ^ _ ^4 

2x V 2 



(47r) 3 x 2 v/i -Aa 
As a; is large, the integral can be performed and one has 



(A12) 



-2x 



t / n • 2/i 3 ^ /, 3 \ e 



(A13) 
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The potentials are given by Eq. |A1|: 



v c (x) - (my£i 

^ W ~ \2m) x^ 



48 a/? m 2 



V, 



LS, 



X 



\2mJ 



(47r) 3 fl 
,-2x 



1 + 



2xJ X^yfx 



16x 



1 | 15 21 \ 

16x 64x 2 / ' 



(A14) 
(A15) 



The numerical value of the terms within square brackets is 275.299 for the central poten- 
tial and 5.88607 for the spin-orbit potential, and they should be compared with those given 
in Eqs. fl5|) and (0), respectively. This result motivated the structure of the expressions 
used in the parametrization of the other components of the potential. 
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FIGURES 

FIG. 1. Diagrams contributing to the minimal chiral background for the 7T7tE-NNP. The vertices 
with one and two pions may be given by either PS and S or PV and V couplings [p5| . The last 
diagram represents the subtraction of the iterated OPE from the reducible box diagram. 

FIG. 2. Ratio between the non-OPEP part of the central (T = 1, S = 0) component of the 
dTRS ||] (dashed line), Argonne vl4 Q (dot-dashed line) and Paris || (continuous line) potentials 
divided by <? 4 and the minimal chiral background for the 7T7tE-NNP. 

FIG. 3. Profile functions for the spin-isospin symmetric component of the dTRS Q (dashed 
line), Argonne vl4 [||] (dot-dashed line), Paris || (continuous line), and the chiral background 
(dotted line) potentials. 

FIG. 4. Same as Fig. || for the (T = 0, S = 1) spin-orbit potential. 

FIG. 5. Same as Fig. || for the (T = 1, S = 1) tensor potential. 
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FIGURE 2 



from paper "Chiral background for the two pion exchange nuclear potential: 
a parametrized version". - C.A. da Rocha and M. R. Robilotta. 




FIGURE 3 

from paper "Chiral background for the two pion exchange nuclear potential: 
a parametrized version". - C.A. da Rocha and M. R. Robilotta. 




FIGURE 4 



from paper "Chiral background for the two pion exchange nuclear potential: 
a parametrized version". - C.A. da Rocha and M. R. Robilotta. 




FIGURE 5 

from paper "Chiral background for the two pion exchange nuclear potential: 
a parametrized version". - C.A. da Rocha and M. R. Robilotta. 



